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Abstract. - We study the spatial correlation function of the phase and its derivative, and re- 
lated, fluctuations of topological charge, in two and three dimensional random media described 
by Gaussian statistics. We investigate their dependence on the scattering mean free path. 



Introduction. - In disordered media the scattering mean free path (notation £) denotes 
the characteristic distance that waves need to travel to achieve a random phase [1] . The length 
£ sets the scale in all aspects of multiple scattering. "Mesoscopic" samples have sizes larger than 
£ and smaller than the decoherence length, beyond which the phase is irreversibly destroyed, 
by either coupling to environment (quantum) or due to background noise (classical) . A third 
fundamental length, the transport mean free path (notation £*), determines the randomization 
of the direction of propagation. It is closely related to £ but the physics is subtly different: 
whereas £ is determined by superposition, £* is also sensitive to interference. 

Despite its conceptual importance, £ is rather unaccessible. Theoreticians have found it 
hard to calculate for dense media because infinitely many "Feynman diagrams" contribute [1], 
and choices have usually been guided by effective medium arguments [2] . In experiments it is 
mostly £* that appears in observables associated with energy transport [3]. To get £ directly 
we must measure the 'coherent' field fi>), either for different sample thicknesses or for different 
times. This requires a coherent source, phase sensitive detection of an exponentially small 
signal, and a good ensemble average. Today, these requirements can be met for ultrasound [4] 
and electromagnetic waves [5] . A final complication is that coherent beam experiments always 
measure the combination of phase randomization and absorption, quantified by the extinction 
length £ e . Here, we investigate the possibility to retrieve the scattering mean free path of a 
random medium from the spatial phase statistics of the scattered waves, rather than from the 
coherent field. The phase of a diffuse field can be measured using a "passive" source, whose 
coherence, position and magnitude are of no relevance, provided it delivers diffuse energy well 
in excess of the background noise. 
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First we introduce the field correlation function in random media and discuss the role of 
small absorption. Secondly we discuss the basic results related to phase statistics and give a 
rapid theoretical support. We finally discuss the relation to topological charge, whose variance 
has been studied extensively in literature. 

Spatial Field Correlation. - In a multiply scattering medium far from the localization 
threshold, it is widely accepted that the complex field obeys circular Gaussian statistics [7]. 
Since the average field (^(r, t)) vanishes if the source is far away, the only parameter governing 
the wave statistics is the two-point correlation function (\&(r — |x, t — ir)\f , (r + |x, t+^r)) [6]. 
We consider a heterogeneous medium in either 2D or 3D, in which all waves suffer from the 
same small absorption rate l/r in energy (absorption length i a — cr ).We assume the medium 
to be large enough to ignore the boundaries. 

We consider a wave field , 0s(r, t) released by a source with spectral density S(oj). The 
subscript B refers to a band-filtering in some frequency band B in which diffusion constant 
D and absorption time r a are constant (Bt > 1 is imposed to resolve the dynamics). Its 
(ensemble-averaged) spatial correlation at any lapse time t in the coda takes the following 
diffuse form [8], 

/ s t _\ 

(* B (r ix,t)* fl (r + ix,<)*^ ~ CXP y^ /2 J b ^S(u) Im G(x,c,C 1 - O 

This formula shows that, once absorption has been acknowledged by the factor exp(— i/r a ), 
the field correlation function is proportional to (the imaginary part of) the absorption-free 
Green function. This is best understood in the frequency domain, where Eq. gives rise to 
the following product of exponentials 

exp[i(w + i/2r a )(t — r/2)] x exp[i(oj + i/2r a )(t + t/2)]* = exp(— iojr) exp(— t/T a ), 

from which infer that the absorption only appears as a function of the lapse time t, and not 
of the local correlation time r. A first observation is thus that in the diffuse time-tail, we can 
consider the (normalized) spatial correlation of the field to be free of absorption. 

In this work we shall ignore the effects of source spectrum and finite bandwidth, which 
can both straightforwardly be included into Gaussian statistics [18]. If we denote by C{x) 
the spatial field correlation function at frequency u [9], in the vicinity of some point r, deep 
in the coda at time t, and normalized to C(0) = 1, we find it to be solely dependent on the 
local Green's function, a central issue in recent attempts to image passively [10]. Contrary to 
chaotic media, the C(x) of random media always decays exponentially due to dephasing from 
scattering. In 3D is C(x) — sinc(fcx) exp(— x/2£), whereas in 2D C{x) = Jo(kx)exp(—x/2£) 
[11], with k = 2ir/\ the wave number. The damped oscillations on the scale of the wavelength 
A originate from a superposition of plane waves incident with arbitrary directions but with 
equal amplitude. It prevents us from seeing the genuine dephasing that occurs on the much 
longer scale of the scattering mean free path. 

Phase correlation. - Rather than the field correlation, we propose here to consider the 
spatial phase correlation Sb(x, t) = ($(r — |x,t)$(r + ix,i)). The phase $(r,i) is defined 
as the complex phase of the wave function = Aexp(i$). At any point x the phase is 
a flat random number between — ir and ir. Upon moving to another point x' it will exhibit 
discontinuous jumps of 2tt between ±7r. An occasional jump of a tt could occur when a 
singularity with either positive or negative topological charge [12] is crossed, but such event 
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Fig. 1 - Spatial cumulative (unwrapped) phase correlation Su(x) versus distance x for a random 
medium in two dimensions. The distance has been measured with respect to the scattering mean free 
path £, and the correlation function has been normalized to k£/n which corresponds to (minus) the 
asymptotic value for large kR. The curves for kl — 10, 100, 200 overlap and follow an exponential 
decay (in red). 

Fig. 2 - Spatial cumulative phase correlation 3u(x) for a random 3D medium, for different values of 
kl. The asymptotic value varies logarithmically with kl. Dashed line is obtained by integrating just 
the asymptotic formula 

Fig. 3 - Cumulative phase variance (&fj(x)} in 2D random media, divided by the optical path kx, 
against distance x, measured in units of the mean free path. Two length scales exist. An initial linear 
rise is visible on the scale of the wavelength. Subsequently, the phase variance reaches the diffuse 
plateau after a few mean free paths. The plateau value depends logarithmically on the value kl. 



has zero probability on a ID line. Let us now "unwrap" the phase, i.e correct by hand for the 
2n discontinuities as was done in the frequency domain [13] and in the time domain [14], to 
get the unwrapped phase or cumulative phase that varies continuously with x, and which 
in principle can take all values [—00,00]. This continuation is not topologically invariant, 
but depends on the exact path chosen to go from r to r + x. Gaussian statistics demonstrate 
that both the distribution P(d&/dx), and the spatial correlation function C$(Aa;) of the phase 
derivative d^/dx are smooth, confirming that phase jumps do not affect phase statistics within 
an infinitesimal increment dx. So the integral of d^/dx along the path should - statistically 
speaking - be equal to the continuous cumulative phase. If we agree to walk straight along 
the z-axis, this results in, 

Eu(x)=J dx' J dx " \^( x ')^ x ")/ = -J dx'x'[C^{x') + C 9 {x-x')). (1) 

The second equality uses that the phase derivative correlation function C$(x' ,x") depends 
only on \x' — x"\ 

The variance of the unwrapped phase has already been discussed in literature in view of its 
close connection with fluctuations in the positions of zeros of the wave function [15, 16, 18, 19]. 
It follows from Eq. -QJ that, 

($l(x))=2[ dx' {x-x')Os,{x'). (2) 
J 

The phase derivative correlation function C<${x) can be calculated from the joint probabil- 
ity function of 4 complex fields at 4 different positions, which involves both C(x), C'(x) and 
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C"(x). The result was earlier reported as an integral [20], but can be put into closed form, 

2 r 

C<3,(x) = - d(f> 
n Jo 

= i(logC)"log(l-C 2 ) (3) 

with H(y) = arctan(^/l + y/y/1 — y)/ yl — y 1 ■ Already for x > A, C is small enough to 
approximate C$(x) by its asymptotic limit, and we find 

L ^ X >X) ^ 2 K } ~ C GJ -\ (l/2x 2 )exp(-x/^) (3D) [A > 

We see that, contrary to the field correlation C, C$ exhibits no oscillations on the scale of 
the wavelength. Equation (J3J shows that the mean free path I assures convergence of Eq. 
in both dimensions, making I the characteristic length scale in spatial phase correlation. In 
fig. 1 we plot as a function of xji for 2D. For all values of kl it follows a universal 

exponential decay towards the asymptotic value —kl/n. In 3D the asymptotic value varies 
only logarithmically with kl in 3D (fig. 2) and two characteristic length scales appear, A and 
I. A measurement of Eu(x) would thus give us direct access to ki and i in both cases. 

The same study has been performed for the phase variance. We infer from Eq. that 
the parameter (&y(x))/x is less sensitive to large values for x than the correlation Eu(x). 
In 2D random media (fig. 3), its plateau value £>$ (the "phase diffusion constant") depends 
logarithmically on kl, although the diffuse plateau is still reached after typically one mean 
free path. In 3D (dashed line in fig. 4) we see that after an approximate linear rise, (&^,(x)) /x 
converges already after a few wavelengths to the asymptote Dq> = 1.922 x k, independent of 
kl. We conclude that the cumulative phase correlation is sensitive to the mean free path in 
both 2D and 3D, unlike the phase variance which depends on I only in 2D. 

Statistics of Topological Charge. - It was realized recently by Wilkinson [16] that the 
spatial phase correlation function is a key element in a old discussion on the statistics of 
topological charge of Gaussian field inside a closed contour. This discussion goes back to 
Halperin [17] and Nye and Berry [6], with significant contributions later by Berry and Dennis 
[18], and Freund and Wilkinson [15]. The zero's of a complex function in d dimensions 
are located on d — 2 "nodal" hypersurfaces, nodal lines in 3D and nodal points in 2D. The 
topological charge Q enclosed by a 2D surface is defined as the number of nodal points weighed 
by their topological sign. This number is determined by the order of the vortex surrounding 
the nodal point [12]. For complex Gaussian random fields topological charges different from 
±1 arc highly improbable [6]. The following relation holds, 

j> dv ■ V$(r) = 2nQ (5) 

with r the line contour enclosing the surface. This relation is reminiscent of the quantized 
rotation of superfluids, and is deeply related to the Rouche theorem in complex analysis. Its 
validity can be understood by applying Stokes' theorem to the function Vlog^(r), whose 
imaginary part has been analytically continued along the contour. 

Equation © facilitates a study of the statistics of the topological charge Q. In the diffuse 
regime the phase gradient in any direction vanishes on average, so < Q >— 0. Following 
Wilkinson [16] we will relate the variance of Q to the correlation of the phase derivative. 
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Fig. 4 - Fluctuations of the topological charge accumulated inside a circle of radius R of a Gaussian 
random field. On the vertical axis is plotted the variance of (Q 2 ) divided by 2nR/X, i.e. the cir- 
cumference in units of the wavelength. Upper two solid lines: 3D random medium with k£ = oo, 10. 
Dashed curve is the cumulative phase variance - divided by 2n, see Eq. @ - obtained by walking 
along a straight line with length 2ttR. The bottom solid curve is calculated for a field correlation 
function C(x) = 2J\{kx)/kx [12]. In all cases the characteristic length scale is the wavelength, and 
not the mean free path 

Fig. 5 - Fluctuations of the topological charge in 2D, accumulated inside a circle of radius R. On 
the vertical axis is plotted the variance of {Q 2 ) divided by the circumference 2tyR/X in units of the 
wavelength. On the horizontal axis the radius of the circle is normalized by the scattering mean free 
path. The characteristic length scale in 2D is the mean free path. The variance of the topological 
charge depends logarithmically on the value hi. 

Fig. 6 - Variance of the smoothed topological charge in 3D (normalized to kl) and in 2D (normalized 
to \ogki). The smoothing is done on Q(R) with the Gaussian exp(— R 2 /R 2 ), as a function of Ro in 
units of the scattering mean free path. For Ro > I it converges to a constant number dependent on 
dimension. 



More precisely, if x R measures the length along a circular contour with radius R we find, 

(Q 2 ( R )) = 77TT2 I de I de'C^(A6) = ^- [ dA0C*(A0) (6) 
{2tt) 2 J J 2ir J a 

The phase derivative correlation C$(A0) = (dg^(9)de^(6')) is, mutadis-mutandis, given by 
Eq. (with (') = d Ag and C = C[2kRsm(A9/2)]). One basic feature of (Q 2 (R)) is already 
known. If we would assume all nodal points n to have random charges Q n = ±1, independent 
of each other, we would find Qn) 2 ) ~ R 2 , i-e. proportional to the surface. Actually, they 
are not independent but tend to be perfectly screened [17,18], which affects the quadratic law. 
We shall define the screening length £ as the typical length beyond which the asymptotic form 
of the charge variance is reached. Wilkinson and Freund report a linear, "diffuse" asymptotic 
form [12], 

B,^oo R TT Jo 1-C(X) 2 1 ' 

We have verified that Eq. © is consistent with this diffusion law, with the same phase diffusion 
constant obtained earlier for the phase accumulated along a straight line (Fig. 3 and 4). Berry 
and Dennis [18] calculate (Q 2 (i?)) from the topological charge correlation function, which is, 
due to screening, sensitive to "edge effects" . When the charge distribution is smoothed with 
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the Gaussian exp(— R 2 /Rq) (with smoothed surface itRq), the quadratic law valid for small 
i?o turns over into the asymptote, 

f°° C(r) 2 
Urn (Q 2 am (Ro))= / dxx [ 
-Ro— oo Jo 1 — C[xy 

provided the integral converges. In figs. 4 and 5 we have calculated (Q 2 ) for a disk in 2D 
and 3D random media. For small kR we see that (Q 2 ) ~ fc 2 i? 2 for both dimensions. In 3D 
(Q 2 ) at large R depends very weakly on the mean free path £ since the integral in Eq. Q 
converges even for £ = oo. For the field correlation C(x) = 2J\{kx)/kx [12] the curve is very 
similar though with a somewhat smaller phase diffusion constant. In both cases this yields 
a screening length proportional to the wavelength. For 2D random media we can see that 
(Q 2 ) depends logarithmically on k£. The topological screening length now depends on the 
mean free path. The variance of smoothed charge (Q 2 m ) however (fig. 6), converges to k£/ir 
in 2D and to 0.54 logk£ in 3D, with a screening length typically equal to the £. Note that 
in 2D chaotic media, for which C(x) = Jo (for), the integrals Q and ((HJ would diverge and 
the topological charge variance would actually increase as (Q 2 ) oc kR\og(kR), i.e. is slightly 
superdiffuse, whereas its smoothed version (Q 2 m ) no longer reaches a constant but becomes 
diffuse: (Q 2 m ) -> kR /(2^). 

Conclusions. - We have shown that in the late coda of waves propagating in 2D and 3D 
random media, the scattering mean free path governs the spatial fluctuations of both the phase 
derivative d^/dx, cumulative phase and the topological charge, measured inside a surface of 
size A. Screening makes topological charge fluctuations grow slower than A. The screening 
length is defined as the typical length beyond which they reach their final behavior. In 2D it 
is proportional to the mean free path, in 3D to the wavelength. If the charge is smoothed with 
a Gaussian the screening length is typically equal to the mean free path, in both dimensions. 
This highlights the subtle role of smoothing screened charge. 

The spatial correlation of d$>/dx, whose measurement requires only four phase-sensitive 
detectors, decays exponentially with the scattering mean free path i as the sole length scale. In 
seismic measurements, often dominated by 2D Rayleigh waves, this may be a novel opportunity 
to measure the scattering mean free path of surface waves. Such measurement would not be 
sensitive to absorption, neither to the precise source location, quite opposed to a measurement 
of the coherent beam. 
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